Unified QCD evolution equations for quarks and gluons 
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Abstract 

Considering the BFKL and DGLAP QCD evolution equations for structure functions, we discuss 
the possibility of unifying them in the whole x and Q 2 range. We emphasize that the main problem 
is related to the constraint of angular ordering of the radiation, and the cancellation of the related 
collinear singularities for inclusive processes. At the leading log 1/x and log Q 2 level, we write 
down a unified system of equations satisfying this cancellation constraint. At low x, it leads to a 
less singular behaviour of the structure functions than the BFKL prediction. 

Resume 

Partant des equations dites BFKL et DGLAP de CDQ pour les fonctions de structure, nous 
discutons la possibilite d'unifier celles-ci dans tout le domaine en x et Q 2 . Nous montrons que 
la principale difnculte vient de l'ordonnancement angulaire des radiations, et de l'annulation des 
singularites collineaires associees dans les processus inclusifs. A l'ordre dominant en log X/x et log 
Q 2 , nous ecrivons un systeme unifie d'equations satisfaisant cette contrainte. Cela conduit a un 
comportement des fonctions de structure a petit x moins singulier que celui predit par BFKL 



Introduction 



k> There exist two different tools in perturbative QCD 
Ih wrien one wants to get predictions for quark and 
C$ gluon structure functions measured in Deep-Inelastic 
Scattering: The first one, valid when Q 2 is large , is a 
resummation to all order of the leading logarithms of Q 2 
(LLQ 2 ), namely terms of the type a n s ln(Q 2 / A.Q CD ) n , 
due to the collinear singularity of the radiative 
corrections. It leads to the so called " Altarelli-Parisi" 
(or DGLAP) evolution equations ^] at leading order. 
The second one, valid when the Bjorken variable x 
is small, is a resummation to all order of the leading 
logarithm of x (LLl/x), namely terms like a"(lnl/x) n , 
due to the infra-red singularities of the radiative soft 

* Talk given in the Proton and Photon structure and diffractive 
interaction session at the Workshop on Deep Inelastic scattering 
and QCD, Paris, April 1995 



part. This tedious calculation was performed twenty 
years ago by L.N.Lipatov and collaborators (BFKL 
equation) and predicts a singular behaviour of the 
proton structure function at small x ||. HERA revived 
this result, since the data on the quark structure 
function inside proton at very small x are in qualitative 
agreement with the BFKL prediction B. 

Since the HERA experiment covers a very large 
range in x and Q 2 , it would be of great interest to have 
a unified system of these two equations. In the past, 
this possibility has already been discussed. It was first 
noticed that a system of equation taking into account 
both LLl/x and LLQ 2 terms can be written down, 
by a precise combination of the corresponding integral 
kernels || . More recently, it has been proved rigorously 
that a unified description of the gluon radiation in 
the whole arrange is possible, due to the property of 
angular ordering p|. In this picture, one can show that 



the dominant contribution due to collinear singularities 
(present in any gluon production amplitude) comes 
from the regions satisfying the following kinematical 
property: 

Q/x > ... > 9i » 9i-i > ... > Ox (1) 

where 6{ ss (qt)i/%i ar e the angles of the emitted 
gluon with respect to the direction of the first emitted 
gluon momentum. The two previous regimes are 
recovered in two different limits: i) at x of the order 
of unity, corresponding to non strongly ordered Xi, one 
recovers the qt ordering and the LLQ 2 resummation . 
ii) at small x, the gluon momentum fractions 
strongly ordered, and the relation ([j]) doesn't imply q t - 
ordering. Thus collinear singularities as well as infrared 
ones can contribute to LLl/x singularities. 

The key point is that for structure functions (this is 
not the case for non-inclusive quantities) the collinear 
singularities cancel in such a way that one recovers the 
BFKL evolution §, [§. We propose here a scheme of 
such unification for both glue and quarks, which fullfill 
this constraint in an effective way B. 

2. Unified QCD evolution equations 

2.1. Equations in Mellin space 

We will restrict ourselves, for the sake of simplicity, 
to the case of a fixed coupling constant as, as for 
the original BFKL derivation. We consider the double 
inverse Mellin transform of the singlet (F s ) and gluon 
(Fq) structure function with respect to Q 2 and x: 



J 2l7T 

x J |e(H w V G(j ; 7 ). (2) 

The DGLAP equations (for fixed as) can be written in 
matrix form for ips and tpQ as follows: 




vg 



, (3) 



where {vg, Vf, 4>gi ^f} are the usual (j-dependcnt) 

DGLAP weights &, and (fs Q are ^ e m ^ial conditions. 

The equation (|3|) must be modified so as to take into 
account the BFKL contribution in the gluon sector, due 
to soft singularities. This dominant contribution can be 
expressed as a singularity in the j — plane at the value 
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where 



X (7) = 2V(1) " m - - 7) ; V>(7) = (5) 

is the eigenvalue-function of the BFKL kernel. Closing 
the path of integration in j around the rightest pole 
given by (0), one gets: 



(?) 
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(6) 



where a saddle point method is used and gives the 
dominant contribution at j c — 1/2, corresponding to 
X (7c)=41n2. 

In order to implement this singular behaviour (JsJ) 
in (|^), we replace the gluonic contribution to the 
anomalous dimension by the following M: 



V G{ 3 ) > V*G{3) = 7 Xil) {VG + 



(7) 



where '3/ is an arbitrary function holomorphic in the j- 
plane near j = 1 and below. Inserted in equation (^|), 
this modification provides a system of unified equations 
mixing the DGLAP and BFKL kernels. Indeed, 
inverting the relation (3), after the replacement vq — > 
v G , one gets 



<Pg 



D(jn) 



i 



4tt 7 u F 




(8) 



with 



D(j,7) = l - 



— (v* 
47T 7 { G 
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(4^) (^G-2n^^). (9) 



Before studying this denominator (^|), let us note 
that in principle <\>p also gets a contribution from 
infrared singularities, such that it should be changed 
into 7 xil) { V G + ^1} — where ^x has the same 
properties as This will be discused elsewhere Q. 

The zeroes of D(j, 7) depend on the region in 
the complex j— plane involved in the inverse Mellin 
transform (|2|) , and thus on the region in x one is looking 
at: 

i) for x of the order of unity, as In 1 /a; <C 1, the 
modification M) has no effect, since the zeroes of D(j, 7) 
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are obtained for small values of 7 (of order as)- In that 
limit, one gets from the definition of x(l) '■ 



Xh) » 1/7 + Otf)\ S g kv g + 0(a J ), 



(10) 



at 



and one recovers the ordinary DGLAP equations 
fixed as)- 

ii) When aglnl/a; = G(l), the singular structure of 
the BFKL kernel drives the relevant domain of the 
integration over 7 in (§) near the "critical" value j c — 
1/2. One recovers the singular behaviour compatible 
with the BFKL calculations. Taking the appropriate 
limit j - 1, a s /(j - 1) = 0(1) : 



D(j,7) oc 1 
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(11) 



£.,2. Constraint on the ip function 



The V'(j) function we have introduced seems to be 
arbitrary when we try to unify DGLAP and BFKL 
equations. In fact, because of the precise cancellation 
of collinear singularities that we emphasized in the 
introduction, the collinear singularities arising from 
quark-loop contribution should cancel at small-x. 
Indeed, angular ordering for radiated quarks (Q) can in 
principle modify the LLl/x singularity. Despite the fact 
that this cancellation is well established for gluons, this 
is not proved when "finite parts" are included [fiol - 

We thus impose this cancellation around j = ji, 
namely, considering D(j, 7) at first order in as (a 
complete discussion will be presented elsewhere ||): 



vo{jh) + v F {j L ) 



vf{Jl) 
aN c log 2 



(12) 



An equivalent constraint is to impose that the 
conformal properties of the BFKL kernel Jll| would 
be preserved, namely that the critical conformal weight 
should stay at 7 = 1/2 ||. This is a statement at LLl/x 
level. 

The system of equations ( |l2| ) leads to a decrease of 
the effective dominant singularity with respect to the 
BFKL one ||. This is in agreement with HERA data 

i- 

Note that the energy-momentum conservation can 
be implemented in the system of equations (|^) in a 
consistant way J8|. It imposes that ^(2) = vq(2) = 2. 
This constraint is not related to the previous one at j\ 
and is not responsible in our model for the decrease of 
the j— plane singularity (in contrary to the discussion of 
Ref. @). 



3. Conclusion 

The unification of the evolution equations for structure 
functions is possible, combining the leading-logarithmic 
contributions in both x and Q 2 variables at fixed ag. 
The constraint due to the mixing of LLl/x and LLQ 2 
leads to a shift down of the LLl/x prediction. 

It would be interesting to know the phenomenolog- 
ical consequences of our system of unified equations. 
Work is in progress in order to take also into account 
the running of as and the next-leading order terms. 
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